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Structure Functions in the Nucleon Resonance Region 
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We report on a detailed study of longitudinal strength in 
the nucleon resonance region, presenting new results from in- 
clusive electron-proton cross sections measured at Jefferson 
Lab Hall C in the four-momentum transfer range 0.2 <Q 2 < 
5.5 GeV 2 . The data have been used to accurately perform 
over 170 Rosenbluth-type longitudinal / transverse separa- 
tions. The precision R = ox /ot data are presented here, 
along with the first separate values of the inelastic structure 
functions Fi and Fl in this regime. The resonance longitudi- 
nal component is found to be significant. With the new data, 
quark-hadron duality is observed above Q 2 = 1 GeV 2 in the 
separated structure functions independently. 



The description of hadrons and their excitations in 



terms of elementary quark and gluon constituents is one 
of the fundamental challenges in physics today. Consider- 
able information on nucleon structure has been extracted 
over the past few decades from separations of inclusive 
lcpton-nucleon cross sections into longitudinal and trans- 
verse structure functions. The original experimental ob- 
servation [1] of the smallness of the ratio R = <jl/<?Ti 
the ratio of the contributions to the measured cross sec- 
tion from longitudinally and transversely polarized vir- 
tual photon scattering, respectively, as measured in deep 
inelastic scattering (DIS) provided the first evidence of 
the fundamental spin-1/2 nature of the partons. Since 
that time, separated structure functions have been mea- 
sured in DIS over a wide range of four momentum trans- 
fer, Q 2 , and Bjorken scaling variable x — Q 2 /2Mv, where 
v = E — E' is the electron energy transfer, and M is the 



1 



proton mass. 

The quantity R is expressed in terms of the fundamen- 
tal nueleon structure functions F\ (purely transverse), 
Fl (purely longitudinal) , and Fi (combined longitudinal 
and transverse) as follows: 
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Precision measurements of R are necessary for several 
fundamental measurements. Extractions of the struc- 
ture function F2, or of the purely longitudinal or trans- 
verse structure functions, Fl and Fi, from cross section 
measurements depend on assumptions for R. The un- 
certainties introduced by this assumption arc highly e- 
dependent, where e is the relative longitudinal polariza- 
tion of the virtual photon in the electron-nucleon scat- 
tering process. Uncertainties in the separation of unpo- 
larizcd structure functions also have a direct impact on 
the extraction of the spin structure functions from spin- 
asymmetry measurements in electron scattering. 

Very few measurements of R have been made in the 
nueleon resonance region. Here, the quantity and preci- 
sion of the existing data (prior to this work) was such 
that it was impossible to study either the mass-squared 
(W 2 ) or Q 2 dependences of the separated longitudinal 
and transverse resonant structure. In a resonance ex- 
citation probed at moderate momentum transfer, large 
values of R and, correspondingly, Fl, are posible, due 
to gluon exchanges between the quarks. These effects, 
as well as the longitudinal character of individual reso- 
nances, are accessible via precision measurements of R. 

In pcrturbative Quantum ChromoDynamics (pQCD), 
Fl k a s , and thus R is expected to decrease logarith- 
mically with increasing Q 2 [2-4]. At low Q 2 , Q~ 2 -type 
dynamic higher twist power corrections, interactions in 
which the struck quark exchanges a gluon with one of the 
spectator quarks in the scattering process, are expected 
to be a significant contribution to R. Additionally, kine- 
matic higher twist effects originating from the binding of 
the quarks in the nueleon (target mass effects), are ex- 
pected to be large. It has been previously reported [6] 
that R measured at intermediate Q 2 in the DIS regime 
[1,5-7] is higher than next-to-leading-order pQCD pre- 
dictions, even with the inclusion of target mass correc- 
tions. This enhanced strength in R relative to pQCD was 
argued to be an observation of non-perturbativc higher 
twist effects [6,7]. 

In contrast, recent quantitative studies [8] of quark- 
hadron duality in the structure function F 2 suggest that, 
even in the resonance region, non-perturbative dynamic 
higher twist effects tend to be small for Q 2 as low as 
1 GeV 2 when the structure function is averaged over 
any of the prominent resonance regions. That is, even 
though the function exhibits resonance (hadronic) struc- 
ture, it tends to average to a global scaling curve consis- 



tent with expectations from DIS data and perturbative 
QCD (hence the term quark- hadron duality). 

The results presented here represent the first detailed 
study of longitudinal strength in the full nueleon reso- 
nance region, to investigate nueleon resonance structure, 
nueleon structure function behavior, and quark-hadron 
duality. 

The experiment (E94-110) ran in the summer of 1999 
in Hall C at Jefferson Lab. An electron beam with a near 
constant current of 60 /iA was provided by the CEBAF 
accelerator with seven different energies ranging from 1.2 
GeV to 5.5 GeV. Incident electrons were scattered from 
a 4-cm-long liquid hydrogen target and detected in the 
High Momentum Spectrometer (HMS), over a scattered 
electron angle, 9, range 12.9 < 6 < 79.9. To account 
for backgrounds from 7r° production and its subsequent 
decay into two photons and then electron-positron pairs, 
positrons were measured in the Short Orbit Spectrometer 
(SOS) and also intermittently in the HMS. Other back- 
grounds included electron scattering from the aluminum 
walls of the cryogenic target cell, as well as electropro- 
duced negatively charged pions. Events from the former 
were subtracted by performing substitute empty target 
runs, while events from the latter were identified and re- 
moved by use of both a gas Cherenkov counter and an 
electromagnetic calorimeter. For more details regarding 
the analysis and the standard Hall C apparatus employed 
in this experiment, see Ref. [9]. Additional discussion 
may also be found in Reference [10], a recent analysis of 
clastic scattering data from this experiment. 
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FIG. 1. Representative Rosenbluth plots for the kinematics 
indicated, as described in the text. 

The inclusive double differential cross section for each 
energy and angle bin within the spectrometer acceptance 
was determined from 

d(7 Y C orr 
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dQ.dE' LAOAE' ' 
where AO (AE ) is the bin width in solid angle (scat 



2 



tered energy), L is the total integrated luminosity, and 
Y C orr is the measured electron yield after correcting for 
detector inefficiencies, background events, and radiative 
corrections. The latter include the bremsstrahlung, ver- 
tex corrections, and loop diagrams standard to electron 
scattering experiments. No corrections were made for 
higher order radiative processes involving two photon ex- 
change, since there exists no reliable inelastic calculation 
for such effects. However, the large kinematic coverage 
in e,W 2 , and Q 2 of this experiment demonstrated min- 
imal (less than 2% [11]) non-linear e dependence of the 
reduced cross section. 

The width of the energy bins was varied with the kine- 
matics such that the corresponding bin width in W 2 
was held constant at 0.05 GeV 2 , where W 2 = M 2 + 
2Mu — Q 2 is the reconstructed mass-squared of the re- 
coiling hadronic system. For each energy bin, a weighted 
average cross section over 9 within the spectrometer ac- 
ceptance was obtained after using a model to correct for 
the angular variation of the cross section from the cen- 
tral angle value. In order to minimize dependence on the 
model used to compute both this and the radiative cor- 
rection, the following iterative procedure was employed: 
a model was used to compute the corrections; the data 
thus obtained were fit to obtain a new model; and this 
resultant new model was then employed to recompute 
the original corrections. These steps were repeated until 
the fitting yielded no further changes, ft was also veri- 
fied that the final cross sections were independent of the 
starting model within 0.6%. A positive byproduct of this 
approach is the availability of a new resonance region fit 
which describes the data here presented [12]. 

Typical cross section statistical uncertainties per en- 
ergy bin were less than 1% with systematic errors, un- 
correlated in e, of 1.6%. The total systematic scale uncer- 
tainty in the cross section measurements was 1.9%. The 
full cross section sample consisted of 32 scans across the 
mass-squared range M 2 < W 2 < 4 GeV 2 . Measurements 
at over 1,000 distinct W 2 , Q 2 and e points were obtained, 
allowing for longitudinal / transverse separations to be 
performed at over 170 fixed W 2 , Q 2 values. 

The extractions of purely longitudinal and transverse 
cross sections and structure functions were accomplished 
via the Rosenbluth technique [13], where measurements 
are made over a range in e at fixed x, Q 2 , and the re- 
duced cross section, a r = do /Y = ot + col = <tt(1+ sR) 
is fit linearly with e. Here, Y is the transverse virtual 
photon flux in the electron-nucleon scattering process. 
Both e and Y were calculated from measured kinematic 
variables. The intercept of such a fit gives the transverse 
cross section ot (and therefore the structure function 
Fi(x,Q 2 )), while the slope gives the longitudinal cross 
section ol, from which can be extracted the structure 
functions R(x,Q 2 ) and Fl(x,Q 2 ). Because R is deter- 
mined by the slope of the fit, relative to ot, the uncer- 
tainty in the extracted value of R (and likewise, Fl) is 



dominated by the uncorrelated uncertainties in the cross 
sections versus e. Prior to a separation being performed, 
data covering a Q 2 range of ±0.5 GeV 2 and W 2 range 
of ±0.05 GeV 2 for W 2 < 3.0 GeV 2 and ±0.10 GeV 2 for 
W 2 > 3.0 GeV 2 were brought to a central value using a 
fit. Multiple fits were utilized to assess the uncertainty 
in this step, which was typically less than 3%. Separa- 
tions were not performed if this kinematic correction was 
larger than 60%. Typical example Rosenbluth plots are 
shown Fig. 1. 





1 




|).lfe<Q 2 <1.0GeV 2 






i i i i n|/ i ii i i \|/ i i i i i i i i i i i 


f l 1 


A 


, S n F is 1 .0< Q 2 <2.0 GeV 2 






i , , \|/ , i , , \|/ , i , , , , 


f t 




2.0< O 2 <3.0 GeV 2 




\l/ i i 


i i i \L i i i i \l/ 1 i i i i i i i i i i 






1 

!l 


3.0< Q 2 <5.0 GeV 2 








Ji.i \i/ , i , , \i/ 



1 1.5 2 2.5 3 3.5 4 



W 2 (GeV 2 ) 

FIG. 2. Measurements of R = <Jl/<Jt, as a function of W 2 , 
for the Q 2 ranges indicated. The error bars shown represent 
both the statistical and uncorrelated systematic uncertainties, 
however the former are negligible in comparison to the latter. 

Values obtained for R are plotted versus W 2 in Fig. 2 
for the Q 2 ranges indicated. Because the data were ob- 
tained at fixed spectrometer angles and electron beam en- 
ergies, while varying the central spectrometer momenta, 
each missing mass spectrum covered a range in Q 2 . It 
is clear from the mass enhancements in the data that R 
exhibits resonant structure, and that this variation with 
W 2 can be quite large. This is the first observation of 
such structure, contradicting the typical assumption in 
inclusive analyses that the resonance contribution to R is 
small or negligible [14,15] in the excitation of resonances. 

The almost twenty well-established nucleon resonances 
with masses below 2 GeV give rise to only three distinct 
enhancements in the unseparated inclusive electron scat- 
tering cross section and, of the three, only the first (the 
lowest mass P33(1232), A) state is not a superposition 
of overlapping resonant states. The second enhancement 
region is often referred to as the Sn, as the unseparated 
cross section here is dominated above Q 2 w 1 GeV 2 by 
the ground state 5n(1535) resonance, even though there 
are overlapping candidate states. In R, however, the sit- 
uation appears to be different. There is a prominent peak 
at W 2 s=s 2.2 GeV 2 , somewhat below the Sn dominated 
mass region. 
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While the A resonance region exhibits smaller and 
rather constant values, R is nonetheless not negligible 
here, nor does it appear to diminish over the Q 2 range of 
this experiment. Yet, the spin- flip required for this pos- 
itive parity isospin / = | excitation suggests that it will 
be purely transverse. Recent measurements from Hall B 
at JLab [16] indicate only a small longitudinal resonant 
component. Therefore, it is probable that the R here ob- 
served indicates a substantial background contribution in 
this regime. 

The data also exhibit resonance behavior in the longi- 
tudinal and transverse structure functions separately, as 
shown in Fig. 3. This does not necessarily follow from 
the resonance behavior of the ratio R = olI g t shown 
above, which could be caused by a cross section that is 
resonant in only the transverse channel. Moreover, not 
only do the data demonstrate significant longitudinal res- 
onance structure, but the VF-dependence of F L is larger 
than that of 2xF x above Q 2 « 1 GeV 2 . 

The peak observed in R at W 2 w 2.2 GeV 2 in Fig. 
2 below the Sn is particularly notable in the Fl longi- 
tudinal channel. This mass is close to that of the elu- 
sive Roper resonance, Pn(1440), the electroproduction 
of which is a topic of much interest [17-21]. 
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FIG. 3. The longitudinal nucleon structure function Fl 
(top), and transverse nucleon structure function 2xFi (bot- 
tom), measured in the resonance region (triangles) as a func- 
tion of x, compared with existing DIS measurements from 
SLAC (diamonds). The curves are from Alekhin (dashed), 
and MRST with (solid) and without (dotted) target mass 
effects included, as described in the text. The prominent res- 
onance mass regions are indicated by arrows. The error bars 
shown represent both the statistical and systematic uncer- 
tainties, with the latter being dominant. 

A precise extraction of information on individual res- 
onances, such as transition form factors, from this inclu- 
sive data must involve a detailed fitting study. At lower 
values of Q 2 < 1 GeV 2 , unitary isobar fits like MAID 
[22] give quite definite and accurate predictions based on 
single pion, two pion, eta, and kaon decay channels for 
resonances below W 2 < 4 GeV 2 . At the higher Q 2 values 
of the data presented here, however, multi-pion effects, 
tails of higher mass resonances, and non-resonant com- 
ponents are very significant and therefore such fits tend 
to undercut the data by factors of two or more [23]. 

Also presented with the resonance region data in Fig. 
3 are the predominantly DIS (W 2 > 4 GeV 2 ) data from 
Rosenbluth separations performed at SLAC [7,6]. Where 
overlapping, the two data sets are in agreement, provid- 
ing additional confidence in the achievement of the de- 
manding precision required for these experiments. In all 
cases, there is a smooth transition between the resonance 
and DIS data in both x and Q 2 . 

The curves shown are next-to-next-to leading order, 
from Alekhin [24], including target mass effects accord- 
ing to [26], and from a recent MRST analysis [25] both 
with and without target mass effects according to [27]. 
Both of these are parton distribution function based pa- 
rameterizations of deep inelastic data. However, MRST 
includes data from other sources as well. Alekhin's cal- 
culations are valid only down to Q 2 = 1 GeV 2 . 

It is clear that some prescription for target mass effects 
is required to describe the data. However, for Q 2 > 1 
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GeV 2 , no additional non-perturbative description seems 
necessary to describe the average behavior of the reso- 
nance region. The resonances oscillate around the scaling 
curves. Furthermore, this is true for a range of different 
Q 2 values, indicating that the scaling curve describes as 
well the average Q 2 dependence of the resonance regime. 
This is consistent with quark-hadron duality, and may 
be counted as the first observation of duality in the sep- 
arated transverse and longitudinal structure functions. 

In summary, we have reported results from a first de- 
tailed study of longitudinal strength in the nucleon res- 
onance region. Contrary to most transition form fac- 
tor fit assumptions, the resonant longitudinal compo- 
nent is found to be substantial. Furthermore, the In- 
dependence of the longitudinal structure function is more 
pronounced than the transverse. Significant strength is 
observed between the Su and P33 resonance mass regions 
in the longitudinal channel. Separated measurements of 
the inelastic structure functions F± and Fl are presented, 
allowing for a first observation of quark-hadron duality 
in these distinct functions independently. 
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